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Abstract 

We discuss a system of a nonlinear Kerr-like oscillator externally pumped by ultra- 
short, external, coherent pulses. For such a system, we analyse the application of the 

m 

^^ ' Kullback-Leibler quantum divergence i^[/9||cr] to the detection of quantum chaotic be- 

haviour. Defining linear and nonlinear quantum divergences, and calculating their power 
-'-J spectra, we show that these parameters are more suitable indicators of quantum chaos 

X' 

H ■ than the fidelity commonly discussed in the literature, and are useful for dealing with 

short time series. Moreover, the nonlinear divergence is more sensitive to chaotic bands 
and to boundaries of chaotic regions, compared to its linear counterpart. 

1 Introduction 

Investigations of the problems of chaos in the dynamics of a quantum system are of current 
interest. These problems are particularly important whenever quantum systems are used in 
quantum information theory models, especially in the context of the problem of decoherence 
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processes in the practical realization of quantum computations. Processes of this type induced 
by interactions with the environment or by any imperfections present in a quantum computer 
[H [2], lower the efficiency of entanglement creation between quantum states. Therefore, one of 
the problems concerning chaotic dynamics is to determine whether a system exhibits regular or 
chaotic evolution. For classical systems, the methods for the determination of the boundaries 
between regular and chaotic dynamics have already been developed and are widely used (Lya- 
punov exponents, power spectra or entropy [3J), whereas in quantum dynamics new methods 
must be found. Because of the linearity of the Schrodinger equation, the methods used in clas- 
sical dynamics cannot be used in quantum theory. Therefore, in order to identify features that 
are not related to quantum chaotic behaviour, it is important to consider quantum systems 
which exhibit chaotic behaviour in their classical limit. In particular, certain methods have 
already been mentioned in ^ and are currently discussed in the papers devoted to the prob- 
lems of quantum chaotic evolution. For instance, the theory of random matrices can predict 
the statistical properties of energy level fluctuations [5] and there is an evident correspondence 
between the successive eigenstates of quantum chaotic systems and the eigenstates of random 
matrices [6l [TJ [H HI [10]. There is also a method based on the decay of the fidelity between 
two quantum states (unperturbed and perturbed) [HI [T31 [12] • It has been proven that in a 
region where a quantum system behaves chaotically, fidelity decays exponentially [m |15l [16] . 
Moreover, the entanglement in the evolving system grows when the dynamics of a quantum 
systems reaches the regions of chaotic behaviour [17j. 

In this paper, we propose to apply parameters based on the KuUback-Leibler divergence 
(KLD) as indicators of quantum chaotic behaviour. By discussing a system of a pumped 
nonlinear Kerr-like oscillator, we show that these parameters are not only indicators of quantum 
chaotic dynamics of the system but they also behave chaotically in the classical sense, despite 
their quantum nature. In particular, we discuss the KLD and show that the linear and first 
nonlinear terms in its expansion can clearly show the boundaries between regular and chaotic 
motion in the dynamics of a quantum system. 



2 Quantum divergence 

Even though it is not a true distance measure, the quantum KuUback-Leibler divergence Kkl 
(also termed relative entropy) is widely used in probability theory and information theory for 
comparisons between two probability distributions. Thus, for the continuous random variables 
P and Q described by their densities p{x) and g(a;), the relative entropy is defined as [HI [19]: 



oo 

Kkl[P\\Q]= j p{x)\og^^dx . 



The divergence Kkl is well-defined for continuous variables and it is always non-negative 
{Kkl[P\\Q] > 0), asymmetric {Kkl[P\\Q] t^ Kkl[Q\\P]) and invariant under the transfor- 
mation of parameters. 

In quantum information theory, for the comparison of two density matrices p and a de- 
scribing two quantum states of a system, the quantum counterpart of the Kullback-Leibler 
divergence Kkl is used and defined as follows pO] : 

KKL[p\\(y] = Tr[p{\np-\na)] . (2) 

From the definition (|2]), it follows that the value of quantum divergence is zero only for p = a, 
while positive values measure the extent of difference between the two density matrices and, 
consequently, between the two quantum states. It is also known that this quantum counterpart 
of the Kullback-Leibler divergence, eq.([2]), is singular whenever the reference state is a pure 
state [21]. In order to circumvent this problem, one can introduce the quantum g-divergence 
defined as [21]: 

K,H|a] = D,rr(pV--)U_,i_o, (3) 

where Dg is Jackson's differential operator, which can be expressed as: 

D,m = lisiLim , (4) 

x[q — 1) 



satisfying the identity: 

D, {f{x)g{x)) = (DJix)) g{x) + f{x) {D^g{x)) + x{q - 1) {DJ{x)) {D^g{x)) . (5) 

In the hniit g — )■ 1, the quantum g-divergence tends to the quantum divergence defined in (j2]). 
The quantum q-divergence can be apphed to determine the purity of the states, which is one 
of the most important problems in quantum information theory. 

The importance of using the Kullback-Leibler divergence in information theory results from 
the fact that many other quantities can be interpreted as the results of applying Kkl to specific 
cases. For instance, mutual information can be expressed by means of the Kullback-Leibler 
divergence as given by [22]: 

I{X;Y) = KKLiPiX,Y)\\P{X)PiY)) , (6) 

therefore, it is the divergence between the product of two probability distributions P{X) and 
P{Y) and the joint probability distribution P{X, Y). Moreover, the Kullback-Leibler divergence 
is related to the Shannon entropy via the relation: 

H{X) = \ogN - KKL{PiX)\\Pu{X)) , (7) 

where P{X) is the true probabihty distribution and PuiX) is the uniform probability distribu- 
tion. 

In this paper, we apply the quantum divergence (j2]) as a tool for the analysis of quantum 
dynamics of the system allowing us to determine whether or not the chaotic regions have been 
reached. Several methods are used to determine whether a quantum system exhibits chaotic 
behaviour. One of them is based on the fidelity between two quantum states described by wave- 
functions and has been discussed in the literature pn[T2| [T3]. One of these states is generated by 
the standard mapping procedure, whereas the other evolves under a slightly perturbed evolution 
operator. It has been shown that whenever the system behaves chaotically (in the quantum 
chaotic sense), the fidelity between theses states decays exponentially. Since the quantum 
divergence ([2]) is defined with respect to two density matrices, in this paper we use density 



matrices to describe the evolution of the quantum system. One of the matrices corresponds to 
the evolution of the unperturbed system, whereas the other describes the evolution of a slightly 
perturbed system under study (similarly to fidelity analysis). 

In [121 [13] it has been shown that whenever the dynamics of a system exhibits quantum 
chaotic behaviour, the differences between the analysed states are clearly visible in fidelity evo- 
lution. Therefore, we presume that the analysis of quantum divergence would also demonstrate 
these differences, and consequently, that quantum divergence could also be used to determine 
the boundary between regular and chaotic behaviour in quantum systems. 

3 The model 

For the purpose of analysis of the quantum divergence ([2]) in quantum chaos theory, it would 
be convenient to use a system whose ability to demonstrate quantum chaotic behaviour was 
confirmed by other well-known indicators of quantum chaos. In this paper, we consider a 
system composed of a Kerr-like nonlinear oscillator externally pumped by a series of ultra- 
short coherent pulses. Systems based on Kerr-like nonlinearity are commonly used as models 
for quantum-optical investigations, for a comprehensive review see for instance [231 [21] {and 
references cited therein). Such systems can be a source of various quantum states of the 
electromagnetic field. For instance, Miranowicz et al. [25] were the first to show that the 
discrete superpositions of arbitrary numbers of coherent states (Schrodinger cats or kitten 
states) can be generated by systems involving Kerr-like nonlinearities. These states are one 
of the most commonly discussed states of quantum optics. In general, models with nonlinear 
media have been widely studied in the context of both nonlinear and quantum optics, and have 
been described in various review papers, e.g. [2S1 [27] . 

It is known that the classical counterpart of systems based on Kerr-like nonlinearity can 
demonstrate regular or chaotic behaviour, depending on their characteristic parameters [23 
[23 [30]- Particularly, the influence of the strength of the external pumping of the nonlinear 
oscillator on the nature of the dynamics of the system was discussed in [31], for cases with or 
without damping. 

The system considered here contains a nonlinear Kerr-like oscillator and is externally pumped 



by a series of ultra-short coherent pulses. It is described by the following Hamiltonian: 

H = Hnl + Hk , (8) 

where Hnl represents the nonlinear Kerr-like oscillator, whereas Hk describes the interactions 
with external pumping. In particular, Hnl and Hk can be expressed as (in units oi h = 1): 



Hnl = f(aO a', (9) 



-(aT) a 
H^ = e (at + a) ^ S{t - kT) , (10) 

A;=l 

where x is the nonlinearity constant, e characterizes the strength of the interaction between the 
nonlinear oscillator and the external filed, o) and a are the usual boson creation and annihilation 
operators, respectively. The parameter T denotes the time interval between two subsequent 
external pulses. Here, we assume that T significantly exceeds the inverse of the external field 
frequency. Hence, the series of the ultra-short external pulses can be modelled by a sum of 
Dirac-delta functions. 

For simplicity, our model neglects damping and, in consequence, we can use the wave- 
function approach in order to describe the dynamics of the system. In particular, in order to 
determine the dynamics of the system, we construct the unitary evolution operators on the 
basis of the Hamiltonians flQfTOj) and apply them repeatedly to the initial state of the system. 
Consequently, we obtain a quantum map describing the evolution of our system: 



|^«(^)) = [UnlUkY |^(t = 0)) . (11) 



The operator IJnl describes the unitary evolution of the system between two subsequent ex- 
ternal pulses (corresponding to a time interval T), according to the formula: 

whereas the operator Uk describes the interaction between the Kerr-like oscillator and a single 



ultra-short pulse. This operator can be written as: 

U^ = e-H^'+^) . (13) 

In fact, this method has been applied to the investigations of finite-dimensional quantum state 
generation |32] and to the comparison of the features of quantum and classical dynamics in 
nonlinear kicked systems [31]. At this point one should mention another, very interesting 
method that allows to find the solution of the problems related to the nonlinear oscillators' 
models discussed in [33] . 

In particular, as already shown in [Slj, the classical counterpart of the system considered here 
can exhibit both regular and chaotic dynamics. For fixed values of the nonlinearity parameter 
X and the time interval between two subsequent pulses T and by varying the strength of 
interaction between the system and the external field, one can obtain regular behaviour of the 
system (if the interactions are weak) or chaotic behaviour in the classical sense (if the excitation 
is strong). The values of the parameters leading to a given character of the system's behaviour 
can be read off from the bifurcation diagram. Thus, using the method described in [SI], we 
generated the diagram (Fig.l) showing the nature of the dynamics of the system depending on 
the value of the parameter e. Fig.l shows the values of the real and imaginary parts of the 
complex parameter a which is the classical counterpart of the annihilation operator a. |ap 
is the energy of the system that is the classical counterpart of the considered system. If the 
value of e increases from ~ 0.3 to ~ 0.7, we observe a region of regular dynamics followed by a 
chaotic band, regular window and a region of deep chaos, successively. 

As mentioned in [Ml [35] , the considered system exhibits quantum chaotic behaviour and, 
therefore, it constitutes a good model for testing the usefulness of various parameters for de- 
tecting chaotic regions and regions preceding the quantum chaotic regions. In this paper, we 
intend to find out whether the quantum divergence is sensitive to quantum chaos and whether 
its behaviour changes near the chaos boundary. As mentioned above, apart from the regions 
of clearly regular dynamics and a region of purely chaotic dynamics (a case of "deep chaos"), 
the classical pumped nonlinear oscillator exhibited a small chaotic region (chaotic band) for 
e ~ 0.36. As shown in [35], in this region, the fidelity did not indicate quantum chaos in the 
dynamics of the system. Here, we examine whether or not the quantum divergence is sensitive 



to such chaotic bands. 

4 Analysis of quantum divergence 

Since the definitions (eqns.( l2|3|) ) of the KLD involve density matrices describing the evolution 
of the system, we applied quantum mapping to generate density matrices, subsequently used 
to determine quantum KLD. 

To determine the quantum divergence, we need density matrices for the evolution of un- 
perturbed and perturbed systems. While the unperturbed density matrix can be derived from 
the wave-function generated by means of the quantum mapping (eq.(1TTl)). its perturbed coun- 
terpart can be generated by a similar procedure where we replace the unitary kick operator 
(eq. (IT3|l ) by another operator Uxp, defined by the relation: 

where A^ describes a small perturbation in the laser field strength. In further considerations 
we assume A^ = 0.001. Therefore, the perturbed density matrix can be expressed as: 

Pp{n) = (UmlUkp) \^{t = 0)){^it = 0)\(u+^U+^) . (15) 

The quantum divergence defined by ([2]) is singular when the reference state is a pure quan- 
tum state. We deal with a pumped oscillator without damping and for this reason, instead 
of using the exact definition ([2]) we shall apply the first terms in the series expansion of the 
logarithm function and thus use the relation |36j : 

lnp-lna = y (-1)"+^ ^^^^ - V (-1)"+^ ^-^^^ . (16) 

n=l n=l 

In particular, we shall concentrate on the first and second terms of such expansion in order to 
investigate their usefulness for the analysis of quantum chaotic systems. 

Thus, we first use the terms corresponding to n = 1 (linear terms) and from the definition 



fITBl) it follows that we obtain a linear quantum divergence (LQD) in the form: 

k'kI = Tr [p{p - a)] , (17) 

If we take the first two terms of the expansion (fT6|) we obtain 



kPl = Tr 



p{p'lip-ir-cr + licr-ir 



The divergence thus defined contains both the linear and the first nonlinear term. In further 
considerations, K^"^^ will be referred to as the nonlinear quantum divergence (NQD). 

4.1 Linear quantum divergence 

As proved in [17], the von Neumann entropy (indicating the existence of bipartite entanglement) 
increased rapidly when the initial state of the system exhibited the behaviour characteristic of 
quantum chaos followed by irregular changes (of small amplitude) depending on the number 
of qubits involved in the system. In our considerations we shall utilise the LQD and NQD 
descriptors in order to determine whether the system is chaotic. 

First, we concentrate on the LQD parameter K^j^. For regular motion, the function Kj^j^ 
oscillates regularly (between and 1) in time (described by the number of pulses) (Fig. 2). 
This means that small perturbations in the laser excitation strength (A^) do not influence the 
final state of the pumped oscillator and, in fact, the perturbed and unperturbed states are the 
same quantum states. The same type of oscillations was observed while considering the fidelity 
between the two quantum states discussed in [35] . 

Moreover, we observed the behaviour analogous to that discussed in [35], when the system 
was inside the regular window and the value of e was close to that corresponding to the boundary 
of the region of deep chaos (e = 0.46). Fig. 3 shows that Kj^[^ oscillates and these oscillations 
are modulated by sine-like variations - the nature of these changes in K)^^ differs from that 
depicted in Fig.2. In the long-time limit, Kj^^ ~ 1 was reached after about 6 x 10^ external 
pulses. This proves the importance of analysis in the long-time limit. Moreover, the effect 
of the observed modulations constitutes an example of quantum beats generated in a system 
with Kerr-like nonlinearity. Therefore, when the quantum system was close to the boundary 



of deep chaos, even in the region of regular dynamics, the changes in K]^'^ were no longer of 
the regular oscillatory nature (characteristic of the regular motion of a quantum system) and 
quantum beats appeared instead. 

Furthermore, we analysed the changes in Kj^[^ in the frequency domain. This was par- 
ticularly useful when the final signal was a composition of several oscillations with different 
amplitudes. Whenever the oscillations had a single frequency (as for the regions of regular 
dynamics - Fig.2), plotting the power spectra was unnecessary, as opposed to the regions of 
chaotic dynamics. 

The power spectra of the LQD and NQD can be obtained by calculating the Fourier trans- 
form of Kj^j^ (or Kj^[^) and the square of its absolute value and subsequently normalising it. 

Fig.4 shows the time-dependence of the LQD (Fig.4a) and its power spectrum (Fig. 4b) 
for the system situated deep inside the chaotic region (e = 0.7). Kj^[^ increased with the 
first ~ 2 X 10^ pulses and then exhibited chaotic changes close to unity (Fig. 4 a). The 
time dependence of the linear part of the entropy consisted of several oscillations of different 
frequencies (Fig. 4b). We assume that these frequencies were distributed regularly in the power 
spectrum, since several groups of these frequencies were observed. The combination of these 
oscillations resulted in irregular changes in Kj^[^ (Fig. 4a). Therefore, the exponential decay 
of the fidelity and its further irregular changes of small amplitude [35], as well as the initial 
monotonic increase and further irregular oscillations in the LQD are characteristic features of 
quantum chaos. These irregularities were composed of several groups of oscillations of various 
frequencies (see the power spectrum). 

4.2 Nonlinear quantum divergence 

Since the chaotic behaviour is the domain of nonlinear systems, it can be useful to introduce the 
nonlinear term in the expansion of quantum divergence in our considerations. Fig. 5a shows 

(2) 

the time dependence of the NQD parameter K]^^ for weak external excitations. This situation 

(2) 

corresponds to regular dynamics and we observe slow oscillations of Kj^^ modulated by other 
oscillations of high frequency. In Fig. 5b we see that only a few frequencies were involved in the 
dynamics of the system and we can identify them in the power spectrum. The slowly varying 
changes in Kj^j^ were of the same frequency as those for the LQD Kj^^^ (Fig. 2). The addition 
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of the first nonlinear term introduced fast oscillations. 

For the first chaotic region (chaotic band between two regular regions - e ~ 0.36), the linear 
part of quantum divergence exhibited almost identical behaviour as for e = 0.1, just like the 
fidelity discussed in [5S]. However, the addition of the first nonlinear term to the quantum 
divergence led to visible changes in its behaviour. Although, in general, the character of the 
oscillations was preserved, fast oscillations were more pronounced than for e = 0.1 (Fig. 6a). 
Much more pronounced changes were observed in the power spectrum. In contrast to the cases 
already discussed, instead of single and sparse lines (Fig. 5b), a multi-peak structure spreading 
over the whole spectrum was observed. The time dependence of the nonlinear part of quantum 
divergence was determined by the sum of a large number of oscillations with various frequencies. 
This is typical of chaotic motion. Moreover, by comparing these results with those obtained 
for the LQD K\^i^ (Figs. 2-4), we observed that the results for the NQD differ considerably 
from those for the LQD. Neither the fidelity, nor the LQD were sensitive to this region of 
chaos (which is visible in the bifurcation diagram), suggesting that this region of chaos appears 
in the classical dynamics only. On the other hand, a strictly quantum parameter (such as 
quantum divergence and, consequently, its nonlinear terms) was able to detect the changes in 
the dynamics of the system in this region. Therefore, the NQD seems to be more sensitive to 
chaotic dynamics than the other parameters mentioned here. 

For e = 0.46, corresponding to the regular window in the bifurcation diagram, the LQD 
parameter K\^i^ exhibited changes indicating the similarity with quantum chaotic dynamics in 
the long-time limit; quantum beats appeared (Fig. 3). This was a long-time effect. The NQD, 
even in the short-time limit, exhibited significant changes in the character of the oscillations, 
as compared with those for the smaller values of e. No vanishing oscillations of Kp^i occurred 
(Fig. 7). The changes in the value of Kj^i^ were more complicated and cannot be described 



by any regular functions (Fig. 7a), particularly when looking at the power spectrum of K 



(2) 
KL 



(Fig. 7b). Again, the final shape of the described function is composed of many oscillations 
with different frequencies, however, they are of different character than those described for 
the linear and nonlinear divergence in the case of the first chaotic band. We observed many 
frequencies, however, a certain pattern of grouping was clearly visible. 

For large values of e (deep chaos), the nature of the time dependence of K]^'^ was completely 



11 



different - no regularities were observed (Fig. 8 a). Moreover, the power spectrum (Fig.Sa) 
was similar to that depicted in Fig.Sa. The frequencies of the oscillations constituting the final 

(2) 

time dependence of Kj^j^ were concentrated around several values, however, we observed many 
more groups of frequencies (frequency bands) and the bands were considerably broader than in 
the spectrum derived from the linear divergence. 

Therefore, the addition of the nonlinear term to the quantum divergence allows us to infer 
additional information about the dynamics of the system. The NQD constructed in this fashion 
was more sensitive to the chaotic evolution (identified in the bifurcation diagrams for the 
classical system), especially in regions where the fidelity and LQD indicated no quantum chaos. 

The fidelity and LQD indicated quantum chaotic behaviour in the long-time limit after 
about 10^ external pulses. If we were to introduce damping processes to our model, this time 
would be too long to observe these features. This does not necessarily imply that quantum 
chaos cannot occur in damped systems. From the analysis of the pumped nonlinear oscillator, 
in which damping processes were included (shown in [3T|), it follows that even after ~ 50 



external pulses, no further changes in the dynamics of the system occurred and the system was 
energetically converged. In this case we would not expect the fidelity to reveal the presence 
quantum chaos. As shown in [35j, if the value of the perturbation A^ is very small, the character 
of fidelity decay cannot be used to determine whether or not the system behaves chaotically. 
Therefore, one must seek other signatures of quantum chaos in the dynamics of the system, 

(2) 

ones that would be sensitive to changes even in such cases. The NLQ parameter Kj^j^ obtained 
by the addition of the nonlinear part to the quantum divergence seems to be appropriate here. 
Since the changes in the NQD and its power spectra were visible even when the fidelity did not 
indicate any reaction (and the chaotic behaviour seemed to occur in the classical sense only), 
we suspect that the NQD can be applied as an indicator of irregular changes in the dynamics 
of a quantum system for shorter times than those discussed for the fidelity-based indicators. 

5 Summary 

We analysed the application of the Kullback-Leibler quantum divergence if[p||o"] to the detec- 
tion of chaotic behaviour in a quantum system. In particular, we discussed a system consisting 
of a nonlinear Kerr-like oscillator externally pumped by a series of ultra-short pulses of coherent 
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field. 

Quantum divergence measures the distance between two probability distributions, in our 
case - the differences between two density matrices, one of which was responsible for the 
evolution of the unperturbed system, whereas the other was defined for a perturbed system. 
We showed that these divergence-based parameters can be used to solve the problems of the 
detection of quantum chaotic behaviours in quantum systems. It is known that the signatures of 
quantum chaos are still being developed, however, we believe that the quantum K-L divergence 
can become one of them. In particular, we concentrated on the dynamics of the system without 
damping and used the first one and the first two terms in the expansion of the K-L divergence in 
order to avoid the singularity that occurs for pure states. In the system under study, the linear 
part Kj^j^ of the expansion of the K-L divergence (the linear entropy-like parameter) clearly 
indicated the chaotic behaviour of the system. However, K\^^ was not sensitive to narrow 
chaotic bands and to the boundary of the deep chaos region. This behaviour resembles that 
discussed in [55] for the fidelity. 

If the parameter involving the nonlinear term of the expansion (what we termed nonlin- 

(2) 

ear quantum divergence, Kkl) is considered, the situation changes considerably. The power 

{'2) 

spectrum of K\.i^ serves as an indicator of chaotic behaviour in the evolution of a system even 
for chaotic bands and when the system is close to the border of deep chaos. This was not 
the case with the linear parameter K\^\^. The changes in the behaviour of K^i^ were visible 
when other parameters {fidelity and K\^j} did not indicate chaotic behaviour, despite the fact 
that from the bifurcation diagram, the considered regions could be identified as chaotic in the 
classical sense. In contrast to fidelity, the parameters described here can be useful whenever 
a short-time analysis is required (for instance, for discussions of damped systems). Therefore, 
we believe that the parameters defined here (and NQD K\^i^ in particular) and their potential 
applications pave the way towards the description of quantum chaotic processes. 
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Figure 1: Bifurcation diagram (3D) for the classical counterpart of the discussed model as a 
function of the strength of external excitation e. The parameters are: T = tt, % = 1. 
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Figure 2: Linear quantum divergence Kj^^ versus the number of the pulses from the external 
field. The parameters are: e = 0.1; A^ = 0.001, T = tt, x = 1. 
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Figure 3: Same as in Fig. 2 but for e = 0.46. 
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Figure 4: Linear quantum divergence Kj^j^ versus the number of the pulses from the external 
field - (a) and the power spectrum of Kj^j^ (from 1500 to 10000 pulses) - (b) for e = 0.7; other 
parameters are the same as in Fig 2. 
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Figure 5: Nonlinear quantum divergence K]^^ versus the number of the pulses from the external 
field - (a) and the power spectrum of Kj^j^ - (b). for e = 0.1; other parameters are the same 
as in Fig 2. 
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Figure 6: Same as in Fig. 5, but for e = 0.36. 
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Figure 7: Same as in Fig. 5, but for e = 0.46. 
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Figure 8: Same as in Fig.5 , but for e = 0.7. 
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